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THE  EVALUATION  OF  BESSEL  AND  NEUMANN  FUNCTIONS  THROUGH  THE  USE 
OF  TAYLOR  SERIES  SOLUTIONS  TO  BESSEL'S  DIFFERENTIAL  EQUATION 

I.  INTRODUCTION 

The  solutions  of  boundary  value  problems  expressed  in  cylindri 
cal  coordinates,  for  example  the  sound  field  resulting  from  a plane 
wave  incident  on  a baffle -transducer , are  often  expanded  in  series 
containing  Bessel  and  Neumann  functions.  The  ability  to  evaluate 
such  solutions  depends  on  the  availability  of  numerical  values 
for  these  functions.  Tables  of  Bessel  and  Neumann  functions  are 
presently  available.  However,  these  tables  are  somewhat  limited  in 
regard  to  the  range  of  order  and  argxment,  and  it  is  not  convenient, 
because  of  computer  storage  limitations,  to  store  function  tables 
in  order  to  evaluate  problems  involving  a large  number  of  Bessel 
and  Neumann  functions  (J^(z)  and  N^(z)  respectively).  Consequently, 
in  the  evaluation  of  problems  requiring  Bessel  and  Neumann  functions 
having  a large  range  of  order  and  argument,  an  analytical  procedure 
which  can  be  used  to  evaluate  the  and  N^(z)  for  such  a range 

of  z and  n is  required. 

The  usual  power  series  representations  of  djj(z)  and  N^(z)  may 
be  used  only  if  z and  n are  small.  If  z is  large,  a limited  number 
of  successive  terms  in  the  power  series  may  be  monotonically  in- 
creasing. In  such  cases,  due  to  the  limited  numerical  significance 
available  in  most  computers,  accurate  results  cannot  be  obtained. 
However,  for  large  orders  and  arguments,  approximate  expressions 
exist  for  evaluating  d^(z)  and  Nj^Cz)*  Such  expressions  are  found 
for  example,  in  Watson.^  Unfortunately,  the  regions  in  which  the 
power  series  and  the  approximations  yield  accurate  results  do  not 
include  the  entire  range  of  order  and  argument  required. 

^ The  purpose  of  this  technical  memorandum  is  to  illustrate  a 
procedure  for  calculating  J,  (z)  and  N^(z)  which wLll  extend  the 

'1^  /ftl 

present  computation  capability  to  a somewhat  larger  range  of  n and 
z.  This  procedure  is  based  on  Taylor  series  solutions 
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to  Bessel's  differential  equation  expanded  about  points  near 
the  argument  desired  for  any  particular  and  N^(z).^ 

While  extending  the  computational  capability  presently 
available,  this  procedure  does  not  permit  calculations  over 
the  entire  range  of  arguments  and  orders  not  covered  by  the 
power  series  or  as)miptotic  expansions.  However,  an  improve- 
ment on  this  procedure  is  currently  being  developed  and  will 
be  described  in  a subsequent  technical  memorandum.  This 
further  development  promises  a capability  for  computing 
and  N^(z)  a continuous  range  of  orders  and  arguments, 
including  those  values  of  n and  z not  covered  by  the  previously 
stated  techniques. 
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II.  TAYLOR  SERIES  SOLUTIONS  TO  BESSEL *S  EQUATION 

A.  The  Method 

The  method  used  in  obtaining  solutions  to  Bessel's  equa- 
tion is  based  on  the  Taylor  series  expansion  of  the  solution  about 

a non- singular  point.  Details  of  this  method  are  available  in 

2 

most  texts  on  ordinary  differential  equations. 

Basically,  the  points  to  be  noted  are  as  follows.  First, 
Bessel's  equation  is  linear  and  homogeneous  with  only  one  singular 
point,  that  being  a regular  singular  point  at  z =0.  Therefore, 
solutions  to  the  equation  have  Taylor  series  expansions  about  any 
point  z = a,  where  a > 0;  and  the  series  are  convergent  in  the 
interval  | z - a | < a.  Second,  since  the  equation  is  of  second 
order  and  is  linear,  there  exist  exactly  two  linearly  independent 
solutions . Therefore  the  new  solutions  will  be  linear  combinations 
of  and  \(z)  • 

B.  The  Taylor  Series  Solutions 

Bessel's  equation  can  be  written  in  the  form 

z^  + (z^-n^)y  = 0.  (1) 

dz^ 

Solutions  to  Taylor  Series  Solutions  of  Bessel's  equation  expanded 
about  z = a may  be  written  as 
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2 

(3^+2ax+a^)  + (x+a)  + (x^+2ax+a^-n^)  y = 0, 


dx 

and  the  solution  becomes 


= ) 


m=0 


A x™. 
m 


(5) 


If  this  series  is  substituted  into  Equation  (4)  the  result  is : 
|^2a^A2+aAj^+(a^-n^)A^^  + |6a^A2+6aA2+(l+a^-n^)Aj^+2aAQ^  x 

+ ^ ^(m+2)  (m+l)a^Ajjj_j_2'*‘(2ni+l)  (ni+l)aA^_l_^+(m^+a^-n^)A^ 

1+A  <5  y 

i-l  m-2J 


m=2 


+ 2aA_ 


® _ n 

X - U. 


(6) 


Equation  (6)  is  an  identity  in  x.  Hence,  the  coefficient  of  each 
power  of  X must  vanish.  Thus 


Ap  = - [aA^  + (a^  - n^)A„] 
2a"^  o 


A^  = " ~ g [6aA2  + (1  + a^  - n^)Aj^  + 2aA^] 

6a 


(7a) 


(7b) 


and,  for  m ^ 2, 


m+2 


(m+2)  (m+l)a* 

+ 2aA  , + A o]  • 
m-1  m-2^ 


[ (2m+l)  (m+l)aA^_^^  +(m^+a^-n^)A^ 


m 


(7c) 


A more  convenient  form  of  equation  (7c)  is 


V4 


(m+4) 


^ |(2iiri-S)  (in+3)aA^3  + ( (iirt-2)®+a®-n® 
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Equations  (7)  define  Ag,  A^,  A^j^  ...  in  terms  of  A^  and  A^.  It  is 

therefore  possible  to  define  two  independent  solutions  to  Equation 
(4)  by  an  arbitrary  choice  of  the  A^  and  A^.  One  such  solution 
is  obtained  by  setting  A^  = 1 and  A^  = 0.  This  solution  is  de- 
noted by  and  is  given  by 


p„(^) 


= I 


B 


m 


(z  - a) 


m 


(8) 


m=0 


The  B are 
m 

B = A for  m = 0,1, 2, 3, 4.  . . . (9) 

mm 

Thus 

P (z)  = 1 . [ ] (z.a)2  + (2.^)3 

^ 2a^  6a^ 

00 

+ I (10) 

m=0 

where 


®o  = 1 » 


Bj^  = 0 , 

/ 2 2v 

T,  (a  -n  ) 

“2 

^ 2a"^ 

B =sija£ 

3 6a3  ■ 


and  for  m 2 0 


B 


‘m+4. 


(m+4)  (m+3)a‘ 


+ [(m+2)^  + a^-n^]  B 


m+2. 


|^(2m+5)  (m+3)a 

■}  ■ 


+ + B^ 


5 


(11a) 

(lib) 

(11c) 

(lid) 


(lie) 
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A second  solution  is  obtained  by  taking  = 0 and  A^  = 1.  This 


solution  is  denoted  by  Q (z)  and  is  given  by 

n 


Qn(^)  = 


C„  (z-a)”  . 


in=0 


As  before 


C = A for  m = 0, 1,2,3,^.  • • . 
mm 


Thus 


.2  _2, 


Q^(z)  = (z-a)  “(2^)  (z-a)^  + ~ 


+ L ^m+4 

m=0 


(z-a) 


2a' 
m+4 


where 


Co  = 0, 


c,  = 1, 


C = - •=— 
^2  2a 


c = [2  - (a^-n^)] 

3 ‘ 6a2  ' 


and  for  m ^ 0 


"m+4 


(m+4) 


i 5 \ (2m+5)  (m+3)a  C 

(nH-3)a^  I 

+ [(m+2)2  + a^-n^l  + 2a  | 


The  general  solution  of  Bessel’s  equation  is  thus 


y = F P„(z)  + G Q„(z) 
where  F and  G are  arbitrary  constants. 


.■^1 


-■ — rr= 


— -V* 


(12) 


(13) 


(14a) 

(14b) 

(14c) 

(I4d) 


(I4e) 


(15). 


! I 

i li 
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C.  Convergence  of  the  Series 

It  is  a well  known  result  of  differential  equation  theory 
that  Taylor  series  solutions  are  uniformly  convergent  in  the  in- 
terval 0 < z < 2a,  where  a singularity  of  the  differential  equation 
exists  at  z = 0,  a is  the  point  of  expansion,  and  no  singularities 
exist  in  the  interval  0 < z < 2a.  Hence,  convergence  of  the  method 
is  assured.  What  is  important  from  the  standpoint  of  computation  is 
the  rate  of  convergence  of  the  series.  It  has  previously  been 
suggested  that  for  an  accurate  evaluation  of  the  series,  the  magni- 
tude of  the  successive  terms  of  the  series  should  be  monotonically 
non- increasing,  that  is,  the  sequences  |Bjjj(z-a)™|  and  ICj^^Cz-a)”*!  (see 
Equations  8 and  12)  should  both  be  monotonically  non- increas ing . 
Equations  lie  and  14e  indicate  that  the  above  terms  are  functions 
of  both  n and  a.  Hence  the  desired  behavior  of  these  terms  is  not 
assured  merely  by  choosing  |z-a|  small.  The  value  of  n must  also  be 
considered.  At  the  present  time  no  rigorous  criterion  has  been  de- 
veloped for  choosing  a as  a function  of  n in  order  to  assure  the 
desirable  convergence  behavior. 
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III.  CALCULATION  OF  Jn(z)  AND  N^(z) 

The  solutions  to  Bessel's  differential  equation  described 
in  the  previous  section  can  be  used  to  compute  the  familiar  solu- 
tions Nn(z).  Specifically,  the  procedure  of  computing 

the  values  of  will  be  shown.  This  will  be  adequate  since 

the  procedure  for  obtaining  the  N^(z)  parallels  the  procedure  for 

The  series  expansion  for  converges  very  rapidly  for 

small  values  of  z.  As  z takes  on  larger  values  convergence  be- 
comes slower.  Quantitatively,  the  maximum  value  of  z for  which 
machine  computation  by  use  of  the  series  is  reasonable  is  that 
value  for  which  the  magnitude  of  successive  terms  in  the  series 
are  monotonically  non-increasing.  The  magnitude  of  the  mth  term 
in  the  series  for  is 


\n  + 2m 
;m  + n) I ’ 


/^\n  + 2m  + 2 

)2l 

(m+1) ; (m+n+1) I 


and  the  magnitude  of  next  term  is 


The  ratio  of  the  second  term  to  the 


previous  one  is 


! m + 1 ) + n + 


. The  largest  value  of  the  ratio  for 


any  given  z and  n occurs  when  m = o,  i.e.. 


• For  a monotonically  non- increasing  series  this 
ratio  must  not  be  greater  than  unity,  so  that 

z < 2 + 1 . ( 16 ) 

Although  the  series  for  (z)  is  not  monotonically  non- increasing 
for  this  range  of  z,  it  is  nearly  enough  so  as  to  allow  its  use  in 
these  computations.  If  a is  needed  for  arguments  larger  than 
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those  shovm  in  Equation  16,  the  series  approach  will  not  yield  use- 
ful results.  The  values  for  J^Cz)  for  z > 2 n/u  + 1 can  be  deter- 
mined using  the  and  functions. 

The  P and  Q f.-ictions  have  been  shown  to  be  linearly  in- 
n ^n 

dependent  solutions  to  Bessel's  equation,  hence  the  Bessel  function 
can  be  written  as  a linear  combination  of  P^  and 


< 

II 

N 

Pn(z)  + B 

Qn(^)  • 

(17) 

From  the  definition  of 

s)  and 

Q^(z)  it  is  seen  that 

a)  = 

1, 

(l8a) 

a)  = 

0, 

(I8b) 

Qn( 

a)  = 

0, 

(18c) 

and 

a)  = 

1. 

(l8d) 

Evaluating  in  Equation  (17)  and  noting  Equation  (18)  it  is 

seen  that 


A=J^(a).  (19a) 

The  value  of  J^(®)  from  Equation  (17)  and  Equation  (18)  shows 

that 


B=j;(a).  (19b) 

Equation  (17)  therefore  may  be  written 

J^iz)  = Jn(a)  P^(z)  + j;(a)  Q^(z).  (20) 

The  values  of  Jjj(a)  and  J^(a)  in  Equation  (20)  can  be  obtained  by 
use  of  the  usual  series  up  to  a=2  n + 1 . Subst5  tution  of  this 
value  of  a into  Equation  (20)  yields 
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J^(z)  = J^(2  + 1)  p^(2)  + j;  (2  + 1)  (^(z).  (21) 

There  exists  a neighborhood  about  z = 2 I n + 1 in  which  the  series 
for  P (z),  Q vz)  and  their  derivatives  will  be  monotonically 
non-increasing.  The  values  of  Jj^(z)  in  the  neighborhood  of 
2N/n  + 1 < z < 2 N/n  + 1 + can  be  obtained  through  use  of  Equation 
(21). 

It  should  be  noted  that  must  satisfy  the  relation 

< 2 + 1 (22) 

since  the  radius  of  convergence  of  the  series  is  2 •Jn+l.  If 
J^(z)  is  needed  for  still  a larger  value  of  z than  2 Vn  + 1 + 

J^(2  v/n  + 1 + e^^)  and  (2  \'n  + 1 + e^^)  can  be  computed  using 

Equation  (21).  Taking  a in  Equation  (20)  to  be  (2  n/ n + 1 + ej^)> 
the  result  is 

Jn(z)  = J^(2  ^nrrr  + e^)  p^(z)  + j;  (2  + 1 + e^)  Q^(z)  (23) 

There  exists  a neighborhood  about  z = 2 n/ n + 1 + in  which  the 
series  for  ^n^^^  will  be  monotonically  non-increasing. 

Thus  J can  be  evaluated  for  all  z up  to  2 n/ n + 1 + + e^.  Note 

that  eg  must  satisfy. 


eg  < 2 N/n  + 1 + (24) 

If  dj^(z)  for  still  a larger  range  of  z is  needed,  this  procedure 
may  be  repeated  until  the  interval  2 vn  + l + e^^  + eg  - - - e^  is 
large  enough  to  cover  the  domain  of  interest. 

This  method  of  calculating  and  has  been  programmed  and 
results  have  been  obtained  for  J^(z),  J^(z),  Jg(z),  Nq(z),  Nj^(z), 
and  Ng(z)  for  3 < z < 10  and  increments  in  z of  Az  = 0.5.  These 
results  were  compared  with  six-place  tabulated  values  found  in 
Reference  (1).  The  calculated  values  agree  to  five  significant 
figures . 


